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CLc ABSTRACT 

I ■ Accepted 12/12/11; Received 22/11/10 ^ost Sunyaev-Zel'dovich (SZ) and X-ray analyses of galaxy clusters try to constrain the 
- cluster total mass (MT(r)) and/or gas mass (Mg(r)) using parameterised models derived from 
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both simulations and imaging observations, and assumptions of spherical symmetry and hy- 
drostatic equilibrium. By numerically exploring the probability distributions of the cluster 
parameters given the simulated interferometric SZ data in the context of Bayesian methods, 
and assuming a /3-model for the electron number density ne('") described by two shape param- 
eters j6 and r^, we investigate the capability of this model and analysis to return the simulated 
cluster input quantities via three parameterisations. In parameterisation 1 we assume that the 
gas temperature is an independent free parameter and assume hydrostatic equilibrium, spher- 
ical geometry and an ideal gas equation of state. We find that parameterisation 1 can hardly 
constrain the cluster parameters and fails to recover the true values of the simulated clus- 
ter. In particular it overestimates MT(r2oo) and rg(r2oo) (A^t('"20o) = (6.43 ± 5.43) x 10'^ Mq 
and Tg(r2oo) - (10.61 + 5.28) keV ) compared to the corresponding values of the simulated 
cluster (Mj(r2Qo) - 5.83 x 10'"* M© and 7'g(r20()) - 5 keV). We then investigate parameterisa- 
tions 11 and 111 in which /g(r2oo) replaces temperature as a main variable; we do this because 
/g may vary significantly less from cluster to cluster than temperature. In parameterisation 
II we relate MT(r2oo) and T„ assuming hydrostatic equilibrium. We find that parameterisa- 
tion II can constrain the cluster physical parameters but the temperature estimate is biased 
low ( MT(r2()o) = (6.8 ± 2.1) x 10'"* Mq and rg(r2()o) = (3.0 ± 1.2)keV). In parameterisa- 
tion III, the virial theorem (plus the assumption that all the kinetic energy of the cluster is 
the internal energy of the gas) replaces the hydrostatic equilibrium assumption because we 
consider it more robust both in theory and in practice. We find that parameterisation III re- 
sults in unbiased estimates of the cluster properties {Mj(r2()Q) = (4.68 + 1.56) x 10'** M© 
and Tg(r2Qo) = (4.3 ± 0.9) keV). We generate a second simulated cluster using a generalised 
NFW (GNFW) pressure profile and analyse it with an entropy based model to take into ac- 
count the temperature gradient in our analysis and improve the cluster gas density distribu- 
tion. This model also constrains the cluster physical parameters and the results show a ra- 
dial decline in the gas temperature as expected. The mean cluster total mass estimates are 
also within Icr from the simulated cluster true values: MT(r2oo) - (5.9 + 3.4) x IO'^^Mq and 
Tg(r2Qo) - (7.4 ± 2.6) keV using parameterisation II and MT(r2oo) = (8.0 + 5.6) x IO'^^Mq 
and rg(r2()o) = (5.98 ± 2.43) keV using parameterisation III. However, we find that for at least 
interferometric SZ analysis in practice at the present time, there is no differences in the AMI 
visibilities between the two models. This may of course change as the instruments improve. 

Key words: galaxies: clusters: general - cosmic microwave background - cosmology: ob- 
servations - methods: data analysis 
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1 INTRODUCTION 

Clusters of galaxies contain large reservoirs of hot, ionized gas. 
This plasma, although invisible in the optical waveband, can be ob- 
served in both X-ray and microwave bands of the electromagnetic 
spectrum through thermal Bremsstrahlung radiation and its scat- 
tering of the cosmic microwave background (CMB) respectively. 
This inverse Compton scattering results in a decrement in the in- 
tensity of CMB photons in the direction of the cluster at frequen- 
cies < 218 GHz, and is known as the Sunyaev-Zel'dovich (SZ) 
effect (Sunyaev & Zeldovich 1970; Birkinshaw 1999; Carlstrom, 
Holder, & Reese 2002). To describe the full spectral behaviour of 
the SZ effect, one needs to consider three main components. These 
include the thermal SZ effect caused by thermal (random) motion 
of scattering electrons, including thermal weakly relativistic elec- 
trons, the kinematic SZ effect caused by peculiar velocity of the 
cluster with respect to Hubble flow, and relativistic effects caused 
by presence of the energetic nonthermal electrons in the hot plasma 
of the cluster that are responsible for synchrotron emission of radio 
halos or relics. However, since the last two processes have signifi- 
cantly smaller effects on the overall spectral distortion at cm wave- 
lengths, we only consider the thermal SZ effect in this paper. More- 
over, we ignore the effects of weakly relativistic thermal electrons, 
which are negligible at cm wavelengths. 

A main science driver for studying clusters through their ther- 
mal SZ signal arises from the fact that SZ surface brightness is inde- 
pendent of redshift. This provides us with a powerful opportunity to 
study galaxy clusters out to high redshift. However, estimating the 
physical properties of the clusters depends strongly on the model 
assumptions. In this paper we aim to show how employing dif- 
ferent parameterisations for a cluster model affects the constraints 
on cluster properties. These tasks are conveniently carried out 
through Bayesian inference using a highly efficient parameter space 
sampUng method: nested sampling (Skilling 2004). This sampling 
method is employed using the package Multinest (Feroz & Hobson 
2008; Feroz, Hobson & Bridges 2009). Multinest explores the high 
dimensional parameter space and calculates both the probability 
distribution of cluster parameters and the Bayesian evidence. This 
algorithm is employed to analyse real multi-frequency SZ obser- 
vations made by the Arcminute Microkelvin Imager (AMI), (AMI 
Consortium: Zwart et al. 2008). 

The rest of the paper is organised as follows. In Section 2, we 
describe the AMI telescope. In Section 3, we discuss Bayesian in- 
ference. Section 4 gives details of how we model interferometric 
SZ data. In Section 5, we describe the modelling of the SZ signal 
using both the isothermal /J-model and an "entropy"- GNFW pres- 
sure model. Section 6 outlines the assumptions needed to estimate 
cluster physical parameters and describes how different parame- 
terisations introduce different constraints and biases in the result- 
ing marginalised posterior probability distributions. In Section 7, 
we describe how to generate a simulated SZ cluster in a consistent 
manner for both models, and in Section 8, we present our results. 
Finally, Section 9 simmiarises our conclusions. 



2 THE ARCMINUTE MICROKELVIN IMAGER (AMI) 

AMI comprises two arrays: the Small Array (SA) and the Large Ar- 
ray (LA) located at the Mullard Radio Astronomy Observatory near 
Cambridge. The SA consists of ten 3.7-m diameter equatorially- 
mounted antennas surrounded by an alimrinium groundshield to 
suppress groimd-based interference and to ensure that the sidelobes 



Table 1. AMI teclinical summary. 





SA 


LA 


Antenna Diameter 


3.7 m 


12.8 m 


Number of Antenna.s 


10 


8 


Baseline Lengths (current) 


5-20 m 


18-110 m 


Primary Beam at 15.7 GHz 


2o;i 


5^5 


Synthesized Beam 


« 3' 


« 30" 


Hux Sensitivity 


30 mjy s-i/2 


3 mJy s-'/2 


Observing Frequency 


13.5-18 GHz 


13.5-18 GHz 


Bandwidth 


3.7 GHz 


3.7 GHz 


Number of Channels 


6 


6 


Channel Bandwidth 


0.75 GHz 


0.75 GHz 



from the antennas do not terminate on warm emitting material. The 
LA consists of eight 13-m diameter antennas. A summary of the 
technical details of AMI is given in Table 1 . Further details of the 
instrument are in AMI Consortium: Zwart et al. (2008). 



3 BAYESIAN INFERENCE 

Bayesian inference has been shown to provide an efficient and ro- 
bust approach to parameter estimation in astrophysics and cosmol- 
ogy by offering consistent procedures for the estimation of a set of 
parameters within a model (or hypothesis) H using the data D 
without loss of information. Bayes' theorem states that: 



Pr{&\D, H) : 



PriP\@,H)Pr{&\H) 



(1) 



Pr(D\H) 

where Pr(&\D, H) = P(&) is the posterior probability distri- 
bution of the parameters, Pr(D\&,H) = X(0) is the likeli- 
hood, Pr(&\H) = 7r(0) is the prior probability distribution and 
Pr{D\H) = is the Bayesian evidence. 

Bayesian inference in practice often divides into two parts: pa- 
rameter estimation and model selection. In parameter estimation, 
the normalising evidence factor is usually ignored, since it is inde- 
pendent of the parameters 0, and inferences are obtained by taking 
samples from the unnormalised posterior distributions using sam- 
pling techniques. The posterior distribution can be subsequently 
marginalised over each parameter to give individual parameter con- 
straints. 

In contrast to parameter estimation, for model selection the 
evidence takes the central role and is simply the factor required to 
normalise the posterior over 0: 
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where D is the dimensionaUty of the parameter space. The question 
of model selection between two models Ho and Hi is then decided 

by comparing their respective posterior probabilities, given the ob- 
served data set D, via the model selection ratio 



R. 



Pr(Hi\D) Pr(D\HOPr(Hi) Zi Pr(Hi) 



(3) 



Pr{Ho\D) Pr{D\Ho)Pr{Ho) Zo Pr{Ho) 

where Pr(Hi)/ Pr{H(f) is the a priori probability ratio for the two 
models. It should be noted that the evaluation of the multidime- 
sional integral in the Bayesian evidence is a challenging numerical 
task which can be tackled by using Multinest. This Monte-Carlo 
method is targeted at the efficient calculation of the evidence, but 
also produces posterior inferences as a by-product. This method 
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is also very efficient in sampling from posteriors that may contain 
multiple modes or large (curving) degeneracies. 

4 MODELLING INTERFEROMETRIC SZ DATA 

In the cluster plasma the central optical depth t is typically between 
0.001-0.01 and the temperature T varies from lO'-lO** K. Thus the 
observed SZ surface brightness in the direction of electron reservoir 
may be described as 



61, : 
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7'=7'CMB 

Here By is the blackbody spectrum, Tcmb = 2.73 K (Fixsen 
et al. 1996) is the temperature of the CMB radiation, /(v) = 
(x^j^ - 4^ (1 + 6{x, Te) is the frequency dependence of thermal SZ 



signal. 



, hp is Planck's constant, v is the frequency and ke 



is Boltzmann's constant. 6{x, T^) takes into account the relativistic 
corrections in the study of the thermal SZ effect which is due to the 
presence of thermal weakly relativistic electrons in the ICM and is 
derived by solving the Kompaneets equation up to the higher orders 
(Rephaeli 1995, Itoh et al. 1998, Nozawa et al. 1998, Pointecouteau 
et al. 1998 and Challinor and Lasenby 1998). It should be noted 
that at 15 GHz (AMI observing frequency) x = 0.3 and therefore 
the relativistic correction, as shown by Rephaeli (1995), is negligi- 
ble for k^T^ < 15 keV. The dimensionless parameter y, known as 
Comptonisation parameter, is the integral of the number of colli- 
sions multipUed by the mean fractional energy change of photons 
per collision, along the line of sight 



^ef^ J-oo 



n,{r)kBT,(r)dl 



Pc(r)dl, 



(5) 



(6) 



where n^ir), P^ir) and are the electron number density, pressure 
and temperature at radius r respectively, ctj is Thomson scattering 
cross-section, ttie is the electron mass, c is the speed of light and dl 
is the line element along the line of sight. It should be noted that in 
equation (6) we have used the ideal gas equation of state. 

An interferometer like AMI operating at a frequency v mea- 
sures samples from the complex visibility plane /v(u). These are 
given by a weighted Fourier transform of the surface brightness 
/v(x), namely 



Ay(x)Iy(x) exp(27r!U ■ x)dx. 



(7) 



where x is the position relative to the phase centre, Ay(x) is the 
(power) primary beam of the antennas at observing frequency v 
(normalised to unity at its peak) and u is the baseline vector in units 
of wavelength. In our model, the measured visibilities are defined 
as 



Vy(U}=Sy(U}+Ny{U), 



(8) 



where the signal component, S y(u), contains the contributions from 
the SZ cluster and identified radio point sources whereas the gen- 
eralised noise part, Ny{u), contains contributions from background 
of unsubtracted radio point sources, primary CMB anisotropics and 
instrumental noise. 

We assume a Gaussian distribution for the generalised noise. 
This component then defines the likelihood function for the data 



1 



1 



(9) 



where is the standard statistic quantifying the misfit between the 
observed data D and the predicted data D''(0): 



(10) 
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where v and v' are charmel frequencies. C is the generalised noise 

covariance matrix 



^ ~ ^v,/ ^vy ^vy ' 

and the normalisation factor is given by 

Zjv = {2nf'''-''^'^\C\^'^, 



(11) 



(12) 



where A',.,, is the total number of visibilities. It should be noted that 
since the main goal of this paper is to demonstrate the effect of 
different parameterisations in modelling the SZ cluster signal, we 
ignore the contributions due to subtracted and unsubtracted radio 
point sources so that the non-Gaussian nature of these sources is 
irrelevant. Moreover, the simulations, used in our analysis do not 
include extragalactic radio sources or diffuse foreground emission 
from the galaxy. The efl^ects of the former have already been ad- 
dressed in Feroz et al. (2009), and here we wish to concentrate 
on the different parameterisation of the cluster. We also note that 
foreground galactic emission is unlikely to be a major contaminant 
since our interferometric observations resolve out large-scale emis- 
sion. 



5 ANALYSING THE SZ SIGNAL: yS-MODEL VERSUS 
GNFW MODEL 

As may be seen from equations (5) and (6), in order to calculate 
the y parameter and therefore to model the SZ signal, we need to 
assume either density and temperature profiles (Feroz et al. 2009; 
AMI Consortium: Zwart et al. 2010; AMI Consortium: Rodriguez- 
Gonzalvez et al. 2011) or a pressure profile (Nagai et al. 2007; 
Mroczkowski et al. 2009; Arnaud et al. 2010; Plagge et al. 2010 and 
Planck Collaboration 20 lid) for the plasma content of the galaxy 
cluster. It is also possible to assume a profile for the gas "entropy" 
and then derive the distribution of gas pressure assuming hydro- 
static equilibrium (Allison et al. 2011). Indeed, in general, one may 
choose to model the SZ signal by assuming parameterised func- 
tional forms for any two linearly independent functions of the ICM 
themiod3mamic quantities. 

Following our previous analysis methodology (Feroz et al. 
2009; AMI Consortium: Zwart et al. 2010; AMI Consortium: 
Rodrfguez-Gonzalvez et al. 2011), we first review the application 
of the isothermal jS-model in modelling the SZ effect and extracting 
the cluster physical parameters demonstrating the impact of dif- 
ferent parameterisations on the inferred cluster properties within 
a model. We then repeat our analysis for the Generalised NFW 
(GNFW) pressure profile, first presented in Nagai et al. (2007), to- 
gether with the entropy profile presented in Alhson et al. (2011) 
to model the SZ effect and derive the cluster physical parameters. 
This approach has potential advantages. It not only removes the as- 
sumption of isothermality but also leads to a density profile that is 
more consistent with the results of the both numerical analysis of 
hydrod3mamical simulations (Voit et al. 2003; Nagai et al. 2006; 
Kravtsov 2006; Hallman et al. 2007) and deep X-ray observations 
of galaxy clusters (Pratt & Arnaud 2002 ; Vikhlinin et al. 2006). 
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5.1 Isothermal y6-model 

This model assumes a yS-profile for electron number density (Cav- 
aliere and Fusco- Femiano 1976, 1978) and a constant temperature 
throughout the cluster 



Tc{r) = 



«e(0) 



\W/2 ' 



hi) 

Tj,(r) = constant. 



(13) 



Here ne(0) is the central electron number density, Te is the elec- 
tron temperature, which is assvmied to be the same as the gas tem- 
perature, r„, and fc is the core radius. It should be noted that in 
our model, yS is considered as a free fitting parameter (Plagge et al. 
2010; AMI Consortium: Zwart et al. 2010 ) and is not fixed to for 
example: (fifn) = 2/3 (Sarazin 1988). 

Using this isothermal j8-model, we can then calculate a map 
of the y parameter on the sky along the line of sight by solving the 
integral in equation (5) analytically (Birkinshaw et al. 1999) 



y{s) = yoll + 



(14) 



where yS > 1/3, ^ is the projected distance from the centre of the 
cluster on the sky such that = + and yo is the central Comp- 
tonisation parameter 



^f^a■^k^T^ne{0)rc - |) 



meC^ 
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(15) 



The integral of the y parameter over the solid angle subtended 
by the cluster is denoted by Fsz, and is proportional to the volume 
integral of the gas pressure. It is thus a good estimate for the total 
thermal energy content of the cluster and its mass (see e.g. Bartlett 
& Silk 1994). Thus determining the normalisation and the slope 
of Ksz - M relation have been the subject of studies of the SZ ef- 
fect (da Silva et al. 2004; Nagai 2006; Kravtsov 2006; Plagge et al. 
2010; Andersson et al. 2011; Amaud et al 2010; Planck Collabo- 
ration 201 ld,e,f,g,h). In particular, Andersson et al. (201 1) investi- 
gated the Ksz - I'x scaling relation within a sample of 15 clusters 
observed by South Pole Telescope (SPT), Chandra and XMM New- 
ton and found a slope of close to unity (0.96 + 0.18). Similar stud- 
ies were carried out by Planck Collaboration (Planck Collaboration 
201 Ig) using a sample of 62 nearby (z < 0.5) clusters observed by 
both Planck and XMM-Newton satellites. The results are consis- 
tent with predictions from X-ray studies (Amaud et al. 2010) and 
the ones presented in Andersson et al. (2011). These studies at low 
redshifts where the data are available from both X-ray and SZ ob- 
servations of galaxy clusters are crucial to calibrate the Ksz - M 
relation and such a relation can then be scaled and used to deter- 
mine masses of SZ selected clusters at high redshifts in order to 
constrain cosmology. 

We calculate the Y^z parameter for the isothermal yS-model 
in both cylindrical and spherical geometries. Assimiing azimuthal 
symmetry, Fcyi reads 



YcyiiR) = dl P^{r)2nsAs 

fn^c J-oo Jo 
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(16) 
(17) 
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where R is the projected radius of the cluster on the sky. 

The integrated y parameter in the case of assuming spherical 
geometry Fsph, is given by integrating the plasma pressure within a 
spherical volume of radius r 



nph(r) 



m^c^ Jo 
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Jo 



(18) 



(19) 



It should be noted that there is an analytical solution for the above 
integral provided that the upper limit is infinity and p > I. However, 
since we study the cluster to a finite extent and /3 varies over a wide 
range including yS < 1, we calculate Fsph numerically. 

5.2 GNFW Pressure Proffle 

As the SZ surface brightness is proportional to the line of sight inte- 
gral of the electron pressure, assuming a pressure profile for the hot 
plasma within the cluster to model the SZ effect seems a reasonable 
choice. In this context, Nagai et al. (2007) analysed the pressure 
profiles of a series of simulated clusters (Kravtsov et al. 2005) as 
well as a sample of relaxed real clusters presented in Vikhlinin et al. 
(2005 , 2006). They found that the pressure profiles of all of these 
clusters can be described by a generalisation of the Navarro, Frenk, 
and White (Navarro et al. 1997) (NFW) model used to describe the 
dark matter halos of simulated clusters. The GNFW pressure profile 
(Nagai et al. 2007) is described as 



(tF(-(in 



a\{b-c)/a ' 



(20) 



where Pei is the normalisation coefficient of the pressure profile, 
is the scale radius and the parameters (a, b, c) describe the slopes of 
the pressure profile at r « rp, r > rp and r <i: respectively. We fix 
the values for the slopes to the ones given in Amaud et al. (2010): 
{a,b,c) = (1.0620,5.4807,0.3292). Amaud et al. (2010) derived 
the pressure profiles for the REXCESS cluster sample from XMM- 
Newton observations (Bohringer et al. 2007; Pratt et al. 2010 and 
Amaud et al. 2010) within r^oo- These pressure profiles also match 
(within rsoo) three sets of different numerical simulations (Borgani 
et al. 2004; PiflFaretti & Valdarini 2008; Nagai et al. 2007). They 
thus derived an analytical function, the so- called universal pres- 
sure profile with above mentioned parameters. This profile has been 
successfully tested against SZ data from SPT (Plagge et al. 2010) 
and the Planck survey data (Planck Collaboration 201 Id). 

We calculate the map of the y parameter on the sky along the 
line of sight by solving the integral in equation (6) numerically. 
However, we note that the central Comptonisation parameter yo has 
an analytical solution 



2o-T/'ei'-p 1 n- 



^)r(V) 



mec2 a r(^) 



(21) 



Similarly, to calculate the thermal energy content of the cluster 

within a sphere with finite radius we use equation (18). In this con- 
text, Anaud et al. (2010) have shown that the pressure profile flat- 
tens at the radius of 5r5oo and used this to define the boundary of 
the cluster. One can thus use this radius to define the total volume 
integrated SZ signal. 



WeC^ Jo 



P^{r')A7ir ^dr' 



(22) 
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6 ESTIMATING CLUSTER PHYSICAL PARAMETERS 

To study the physical parameters of the cluster, such as its total 
mass and gas mass, we have to make some assumptions about the 
dynamical state of the cluster. The most widely-used assumptions 
are: that the gas distribution is in hydrostatic equilibrium with the 
cluster total gravitational potential dominated by dark matter, and 
that both dark matter and the plasma are spherically syimnetric and 
have the same centroid. 

The cluster mass Mjirx) is also defined as the total amount of 
matter internal to radius rx within which the mean density of the 
cluster is X times the critical density at the cluster redshift. Math- 
ematically, the assumption of hydrostatic equilibrium applies ev- 
erywhere inside the cluster and relates total cluster mass internal to 
radius r to the gas pressure gradient at that radius and hence to the 
density and temperature gradients respectively 



dlngg(r) dlnrg(r) 



dlnr 



dlnr 



(23) 



where n = 0.6mp (Sarazin 1988) is the mean mass per gas particle, 
wjp is the proton mass and G is the universal gravitational constant. 
Assuming spherical geometry, it is also possible to calculate the gas 
mass and total mass internal to radius rx 



Mj{rx) 



Jr^rx 




ir)dr. 



An 

An . 



(24) 



(25) 



ftntCz) 

universe 



at the 



cluster 



is the critical density 
redshift z and H{z) 



Here 
the 

Ho ^J{aM + nA0(l + zT + nR(l + zf + nK(l + z)^ is the Hubble 
parameter at redshift z, where Hq = 100/jkm s"'Mpc~' is the 
Hubble constant now. CIm measures the present mean mass density 
including baryonic and nonbaryonic dark matter in non-relativistic 
regime, takes into account the present value of the dark energy, 
f^R measures the current energy density in the CMB and the low 
mass neutrinos, fix describes the curvature of the universe and 
g = (1 + 2)3[»'o+"'a] exp [^Y^] is the dark energy equation of state. 

Moreover, it has been long known that the total mass of the 
cluster is strongly correlated with its mean temperature. This arises 
from both X-ray observations of galaxy clusters (Voit & Ponman 
2003) and the fact that the gravitational heating is the dominant 
process in the clusters within the hierarchical structure formation 
scenario (Kaiser 1986 ; Sarazin 2008). 

Assuming virialisation and that all cluster kinetic energy is in 
gas internal energy suggests that T oc M^'', where T is the mean 
gas temperature within the virial radius and M is the cluster to- 
tal mass internal to that radius. However, an extensive range of 
studies based both on observations of galaxy clusters and on nu- 
merical simulations have been carried out aiming to determine the 
proportionality coefficient of such relation (Evrard et al. 1996; Eke 
et al. 1998; Voit 2000; Yoshikawa et al. 2000; Finoguenov et al. 
2001 ; Afshordi & Renyue 2002; Evrard et al. 2002; Sanderson et al. 
2003; Borgani et al. 2004; Voit 2005; Arnaud et al. 2005; Vikhlinin 
et al. 2006; Afshordi et al. 2007; Maughan et al. 2007 and Na- 
gai et al. 2007). Finoguenov et al (2001) studied the observational 
mass-temperature relation of two sets of cluster samples. In their 
first sample they used the assumption of isothermality whereas in 
the second set they knew the temperature gradient of the clusters 
within the sample. In both samples, they found that the discrepancy 
from the self- similarity in the M-T relation is more pronounced 
in the low mass clusters (^B^g < 3.5 keV) as non-gravitational 



processes become more dominant in these clusters. Similar results 
were obtained by Amaud et al. (2005) when they analysed a sam- 
ple of 10 nearby (z < 0.15) relaxed clusters in the temperature 
range 2-9 keV. They showed that the slope of the M-T relation 
for hot clusters is consistent with self-similar expectation while 
for low temperature (low mass) clusters the slope is significantly 
higher. Studies of the observational mass- X-ray luminosity relation 
(Maughan et al. 2007) also show that the scatter in the Lx - A/500 
relation is dominated by cluster cores and is almost insensitive to 
the merger status of the cluster. Theoretical studies based on the 
adiabatic simulations and the hydrodynamical simulations of clus- 
ter formation with gravitational heating only also verify the slope 
of 3/2 in M-T relation (Evrard et al. 1996; Eke et al. 1998; Voit 
2000; Yoshikawa et al. 2000) while nimierical simulations which 
take into account the non-gravitational heating processes and the 
effect of the radiative cooling of the gas (Borgani et al. 2004; Nagai 
et al. 2007) do predict a slightly higher slope. Moreover, almost all 
of the above mentioned studies do agree that the discrepancy in the 
slope of the M-T relation could also be due to the different proce- 
dures used for estimating masses in simulations and observational 
analyses. 

In this paper we therefore decided to follow the approach 
given in Voit (2005). This is based on using the virial theorem to 
relate a collapsing top-hat density perturbation model to a singu- 
lar truncated isothermal sphere. It also takes into account the finite 
boundary pressure and assumes all kinetic energy is internal energy 
of the hot plasma. This gives 



of k^iTgirx) = '^(-j [GM^(rx)H{z)?'\ 



(26) 



It should be noted that above relation assumes that the virialisation 

occurs at rx. 

Based on the above assumptions, one can adopt different pa- 
rameterisations to study physical properties of the cluster within 
a particular model (e.g. using either the assumption of hydrostatic 
equilibrium or the M-T relation). These different approaches shed 
light on the realism of the assumptions made throughout the analy- 
sis and reveal different biases and constraints associated with them. 
In a single- frequency observation that at least partially resolves 
the cluster, the best one can hope to achieve in constraining an em- 
pirical model of the SZ decrement is to estimate the central posi- 
tion of the cluster (the position of the decrement) and two further 
parameters-i.e. shape and scale parameters. The interpretation of 
such constraints does however depend on the particular parrame- 
terisation. 

Hence in the following sections we discuss different pos- 
sible parameterisations within two models; isothermal yS-model 
and "entropy"-GNFW pressure model. In doing so we try to dis- 
entangle the thermal pressure built-in correlation between pairs 
of physical parameters that lead to the SZ effect intensity-i.e. 
(r,,M,),(/,,Mr),etc. 



6.1 Isothermal j8-Model 

Generally, there are two different parameterisations that one could 
use in the analysis of the cluster SZ effect and deriving its phys- 
ical parameters. However, the assumption of isothermality pro- 
vides another form of parameterisation where the gas temperature 
is assimied as an input free parameter along with the assumption 
of hydrostatic equilibriiun (e.g. isothermal >S-model or isothermal 
GNFW model). 

In the following sections we discuss our three different param- 
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eterisations for the isothermal yS-model within our Bayesian frame- 
work. It should be noted that in all of these parameterisations, we 
employ physically-based sampling parameters. Such parameters re- 
veal the structure of degeneracies in the cluster parameter space 
more clearly than parameters that just describe the y-map such as 
angular core radius 6^, shape p and central temperature decrement 
AFo. We also note that throughout our analysis we impose the ad- 
ditional constraint that the cluster has a non-zero r^oa- 



6.1.1 Parameterisation I 

Our sampling parameters for this ease are 0;. = 
(XcjcfcyS, rg,Mg(r2oo),z), where Xc and are cluster pro- 
jected position on the sky, and jS are the parameters defining 
the density profile, Tg is the gas temperature, Mg(r2oo) is the gas 
mass internal to radius r2oo and z is the cluster redshift. It should 
be noted that AMI can typically measure the overdensity radii rsoo 
and ^200 for z > 0.15. However, we choose to work in terms of an 
overdensity radius of r2(x) since the constrains from AMI data on 
the cluster physical parameters are stronger at this radius and this 
radius is approximately the virial radius. We further assume that 
the priors on sampling parameters are separable (Feroz et al. 2009) 
such that 



;r(0,) = n{x,)n{y,)n{r,)n(J3)n{T^)n{M^{r2ao))Mz). 



(27) 



This implies that parameterisation I ignores the known apriori cor- 
relation between the cluster total mass and gas temperature. We 
use Gaussian priors on cluster position parameters, centred on the 
pointing centre and with standard deviation of 1 arcmin. We adopt 
uniform priors on the cluster core radius,/? and the gas temperature. 
As mentioned in Feroz et al. (2009), for SZ pointed observations, 
where we know the cluster redshift from optical studies and pos- 
sibly the gas mass fraction from X-ray studies, we can assume a 
separable prior on the gas mass and redshift, namely, ti{M^,z) cc 
7i{M = Mg/fg,z) = n(M = Mg/fg)n(z), (Feroz et al. 2009), where 
each factor has some simple functional form such that their prod- 
uct gives a reasonable approximation to a known mass function 
e.g. the Press-Schechter (Press & Schechter 1974) mass function. 
We will assume such a form in our analysis where MMg) will be 
taken to be uniform in log in the range (M„),^in = 10'' h^^Mo to 



(M„)„ 



5 X 10 h Mq and the redshift is fixed to the cluster 



redshift. A summary of the priors and their ranges for this parame- 
terisation is presented in Table 2. 

Having established our physical sampling parameters, mod- 
elhng the SZ signal is performed through the calculation of the y 
parameter which requires the knowledge of parameters describing 
the 3-D plasma density and its temperature, namely r^p, ne(0) and 
Tg. We sample from , p and Tg but as shown below, deriving 
ne(0) requires employing the assumptions of hydrostatic equilib- 
rium, isothermality and spherical geometry right from the begin- 
ning of the analysis. 

Substituting the isothermal j8-model into the equation of hy- 
drostatic equilibrium equation (23), we can then relate the MT(r2oo) 
to our model parameters as well as to the temperature 



MT(r2oo) ■■ 



3pr 



3 

200 



«2 1 «2 
c ^ '200 



^tG ' 



(28) 



where we have used gg{r) = Hene{r) with //e = 1.14/Mp (Jones et al. 
1993; Mason & Myers 2000) defined as the mean gas mass per 
electron. By combining equations (25) and (28) at r2oo, vve first 
calculate the overdensity radius of r2oo and since Mg(r2oo) is also 



Table 2. Summary of the priors on tlie sampling parameter set in parame- 
terisation I. Note tliat N{}1, tr) represents a Gaussian probability distribution 
witli mean and standard deviation of cr and U(a, b) represents a uniform 
distribution between a and b. 



Parameter 


Prior 






^(0 , 60)" 




''c 


(7(10 , 1000) ft- 


-'kpc 


ys 


£7(0.34 , 2.5) 




logMg(r2oo) 


U(12 , 14.5) r 






£7(0 , 20)keV 





one of our sampling parameters we can recover the central electron 
number density by rearranging equation (24): 



''200 - 



ne(0) ■■ 



9pk^T, 



4;r/<G(200eeri.(z)) 



Ann, 



r 

Jo 



r'^dr' 



(29) 



(30) 



For cluster physical parameters we use the value of the overdensity 
radius of ^200 (equation 29) to calculate the cluster total mass in- 
ternal to r2oo assuming spherical geometry, (equation 25). The gas 
mass fraction at r2oo is then simply /g(r2oo) = Afg(f'2oo)/Wr(''20o). 
As the central electron number density and plasma temperature are 
assumed to be constants, we can in principle calculate cluster phys- 
ical parameters in any overdensity radius other than r2oo by assum- 
ing that the hydrostatic equilibrium holds everywhere in the cluster. 
In this paper we study the cluster properties at two overdensity radii 
r2oo and rsoo. Extracting cluster physical parameters at rsoo in par- 
ticular enables us to compare our results with the results obtained 
from X-ray analysis of the clusters of galaxies. r^Qo is calculated 
by equating equations (23) and (25) and setting X = 500. Mg(r5oo) 
and Mj-(r5oo) are then derived using equations (24) and (25) respec- 
tively. 



6.1.2 Parameterisations II and III 

Parameterisation I does not take into account the correlation be- 
tween the cluster total mass and its mean gas temperature. How- 
ever, as mentioned earlier, observations of galaxy clusters and the- 
oretical studies have both shown that there is a strong correlation 
between these two cluster parameters. We have already used this 
parameterisation in the analyses of 7 clusters out to the virial ra- 
dius (AMI Consortium: Zwart et al. 2010). We found that using 
this parameterisation along side the assumption of isothermality 
led to strong biases in the estimation of cluster parameters. This 
implies that, in the absence of a measured temperature profile, we 
should eliminate gas temperature from the list of our sampling pa- 
rameters and instead sample from either Mj(j'2fya) or Mg(r2oo)!and 
/g(''2oo)- We choose total mass as a sampUng parameter since this is 
consistent with our cluster detection algorithm and analysis (AMI 
Consortium: Shimwell et al. 2010). This form of parameterisation 
then allows us to calculate gas temperature either by using isother- 
mal hydrostatic equilibrium (parameterisation 11), or virial relation, 
(parameterisation III). 

Our sampling parameters for these two parameterisations are 
0c = {Xc, yc, fcP, MT(r2oo), /g(r2oo), z) which are assumed to be in- 
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Table 3. Summary of the priors on the sampling parameter set in parame- 
terisations II and III. 



Parameter 



Prior 



rc 



logMT(r2oo) 
fgirioo) 



N(0 , 60)" 
UilO , 1000)/i-lkpc 
[7(0.34 , 2.5) 
(7(14 , 15.5)/!"'Mo 
N(O.OM , 0.016) /!-' 



dependent for the same reasons described in previous section such 
that 

;r(0c) = 7:ix,) Ttiyc) 7:ir,) niMjir^)) n{f^{rwo)) 7({z). (31) 

The priors on Xc, yc, ''c. yS and z are the same as for parameterisa- 
tion I. The prior on MT(r2oo) is taken to be uniform in logM in the 
range Mn,},, = lO''* /i^'Mq to M^ax = 5x10'^ /i"'Mo and the prior of 
/g(r2oo) is set to be a Gaussian centred at the WMAP7 best fit value: 
/g = 0.12 with a width of 0.016 (Komatsu et al. 2011; Larson et al. 
201 1 and AMI Consortium: Rodrfguez-Gonzalvez et al. 201 1). A 
summary of the priors and their ranges for these two parameterisa- 
tions are presented in Table 3. We calculate Mg(r2oo) from the defi- 
nition of gas mass fraction and r2oo is determined assuming spher- 
ical geometry, equation (25). Central electron number density is 
then calculated using equation (30). 

For parameterisation II, the gas temperature at rioo is esti- 
mated assuming hydrostatic equilibrium using equation (28) and 
is assumed to be constant throughout the cluster 



^Brg(r2oo) 



/jGMT(r2oo)(>foo + rg) 

3/5''200 

(4;rpG)(200e„u(z))(4o + rj) 
9fi 



(32) 



where the last form is derived by substituting for MT(r2oo) us- 
ing equation (25). We refer to relation (32) as the HSE mass- 
temperature relation. Similar to parameterisation I, once temper- 
ature and central electron number density are determined, we can 
calculate cluster physical properties at the overdensity radius of rjoo 
using equations (23), (24) and (25) and setting X = 500. 

For parameterisation III, we calculate the mean gas tempera- 
ture within a virial radius of r2oo using the mass-temperature rela- 
tion described in equation (26), which is then assumed to be con- 
stant throughout the cluster 



feBrg(r2oo) = 8.2(^^j5^J ^— J keV. (33) 

We refer to this relation as the virial mass-temperature relation. 
This also implies that virialisation occurs at r2oo- In our analysis, 
we use this relation to determine the mean gas temperature and 
once this is determined we repeat the same procedure carried out 
for parameterisations I and II to obtain cluster physical properties 
at the overdensity radius rsoo- 



6.2 Entropy-GNFW Pressure Model 

As it was mentioned in Section 5.2 the choice of the GNFW pres- 
sure profile to model the SZ signal is reasonable as the SZ surface 
brightness is proportional the line of sight integral of the electron 
pressure. However, in order to link the gravitational potential shape 



to the baryonic physical properties of the ICM, one has to make as- 
sumptions on the radial profile of another thermodynamical quan- 
tity. Among the thermodynamical quantities of the ICM, entropy 
has proved to be an important gas property within the cluster. En- 
tropy is conserved during the adiabatic coUapse of the gas into the 
cluster gravitational potential well, however, it will be affected by 
any non-gravitational processes such as radiative cooling, star for- 
mation, energy feed back from supernovae explosions and Active 
Galactic Nuclei (AGN) activities. It therefore keeps a record of the 
thermodynamic history of the ICM (Voit 2000; 2003, 2005; Pon- 
mann et al. 1999; Pratt & Arnaud 2002; Allison et al. 2011). 

Moreover, for a gravitationally collapsed gas in hydrostatic 
equilibrium , entropy profile is expected to have an approximate 
power law distribution (« r") (Lloyd-Davies et al. 2000; Voit 
2005; Nagai et al. 2007; Pratt 2010). However, there is a large de- 
viation from self-similarity in the entropy radial profile in the inner 
region of the cluster (r < 0. lr2oo) due to the impact of all of the non- 
gravitational mechanisms described above on the thermodynamics 
of the ICM (Finoguenov et al. 2002; Ponman et al. 1999, 2003; 
Lloyd-Davies et al. 2000; Pratt 2010). In the inner region, the re- 
sults of the non-radiative simulations and simulations that take into 
account AGN activities plus preheating models predict a flat core in 
the entropy distribution due to entropy mixing (Wadsley et al. 2008; 
Mitchell et al. 2009). The observed entropy profiles using X-ray 
telescopes also flatten in the iimer regions in general while having 
similar external slopes (Pratt et al. 2006). In the outskirt of the clus- 
ter (out to virial radius and beyond), on the other hand, the results 
of the latest numerical simulations (Nagai 2011; Nagai et al. 2011) 
and observational studies of the clusters using Suzaku and XMM- 
Newton satellites at large radii including A1795 (Bautz et al. 2009), 
PKS 0745-191 (George et al. 2009), A2204 (Reiprich et al. 2009), 
A1413 (Hoshino et al. 2010), A1689 (Kawaharada et al. 2010), 
Virgo cluster (Urban et al. 2011) and Perseus cluster (Simionescu 
et al. 201 1) also show that behaviour of ICM entropy deviates from 
the prediction of a spherically systematic shock heated gas model 
(Tozzi & Norman 2001). According to these studies major sources 
of this deviation may be due to incomplete virialisation, departure 
from hydrostatic equilibrium, gas motion and gas clumping. 

In this context and to derive the cluster physical parameters 
we decided to adopt the entropy profile presented in Allison et al. 
(201 1) which is a yS- model like profile: 



(34) 



where K^ir) is the plasma entropy at radius r, K^i is the normali- 
sation coefficient of the entropy profile, and a are the parame- 
ters defining the shape of the profile at different radii. Assuming 
an entropy profile with above form guarantees the flat shape in the 
inner region and a power law distribution at the larger radii (up to 
rsoo) where ^^e('') <^ r^" with a ~ 0.55. We note that in order to 
take into account the behaviour of entropy at the clusters outskirts 
we need to modify the assumed entropy profile with additional pa- 
rameter and/or component. However, as the studies of this kind in 
tmderstanding the physics of the cluster outskirts and accurate mea- 
surements of the ICM profiles in the cluster outer regions are still 
ongoing, we do not study a modified form of our assumed entropy 
profile here. We of course aim to consider a more general form in 
our future analyses. 

As for using the GNFW profile, one has indeed to make an 
assumption on either the density, temperature or the entropy pro- 
file shape in order to link the gravitational potential shape to the 
baryonic physical properties of the ICM. In this paper, we decided 
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Table 4. Summary of the priors on the sampUng parameter set m the 
entropy-GNFW pressure model. 



M„(r2oo) 



Parameter 


Prior 




^(0 , 60)" 


fc 


C/(10 , 1000)ft-'kpc 


a 


J/(0.0 , 1.0) 


logMT(r2oo) 


1/(14, 15.5) r'Mo 


/g(f200) 


JV(0.084 , 0.016) 




£7(0.001 , 3)/i-iMpc 



to work with an assumption of entropy profile for all the reasons 
given above. The combination of the GNFW pressure and the " 
a-model" entropy profiles can then fully describe the large-scale 
properties of clusters as they determine the form of the dark matter 
potential well in addition to the structure of the ICM. 

To relate the entropy to the other thermodynamical quantities 
inside the ICM we use the definition of the entropy given in the 
astronomy literature. For an adiabatic monatomic gas, 



= 

Pe = 



(35) 
(36) 



which is related to the true thermodynamic entropy per gas particle 
via 5 = In(A^e) + ■^o where So is a constant (Voit 2005). 

Using equations (20),(34),(35) and (36) we can derive the 3-D 
radial profiles of the electron number density and the temperature. 



ne(r) = nA — 



1 + 1- 



rc 



1 + 1^ 



l(¥) 



rc 



where 



and 



(^ei) 
- p2/5 1 



(37) 



(38) 



(39) 



(40) 



are the normalisation coefficients for the electron number density 
and the temperature profiles respectively. We note that the above 
derived electron number density has components that take into ac- 
count both the fit for the inner slope of the cuspy cluster density 
profiles and the steepening at larger radii (r > rsoo) (Pratt & Ar- 
naud 2002, Vikhlinin et al. 2006). 

As using this model to analyse the cluster SZ signal removes 
the assumption of isothermality, parameterisation I which assumes 
a single core temperature as a free input parameter can not be used 
in the analysis using entropy-GNFW model or any non- isothermal 
model. We therefore study the cluster SZ signal and its physical 
properties using parameterisations II and III. 

Our sampling parameters for this model are ©c = 
{Xc, yc,rc, a, Mr(r2oo), /g(''2oo)> z). A summary of the priors and their 
ranges for the "entropy"-GNFW pressure model is presented in 
Table 4. 

;r(0c) = Ttixo) n(yo) n{,r^) n(a) ;r(Mj-(r2oo)) n{f^{riwi)) i^iz). (41) 

Sampling from MT(r2oo) in both parameterisations leads to the es- 
timation of r2oo assuming spherical geometry for the cluster, equa- 
tion (25). Sampling from MT(r2oo) and /g(r2oo) also allows us to 
calculate Mg(r2oo). «ei is then 



Ann, 



f 

Jo 



1 -I- 



-.(42) 



1 + 



dr' 



In parameterisation II we substitute the electron number density 
and GNFW pressure profiles in the assumption of hydrostatic equi- 
librium (j'--^^ = -^^) at r2oo and derive P^. The normali- 
sation coefficients for the temperature and entropy profiles, k^T^ 
and Kgi, are then calculated by solving the equations (39) and (40) 
simultaneously. 

In parameterisation III we calculate feB7'e(r20()) using virial M- 
T relation, equation(33). kBT^ is then calculated by substituting the 
values derived for r^m and ^B7'e(''2oo) in temperature profile given 
in equation (38). Similarly, the normalisation coefficients for the 
pressure and entropy profiles, and Kd, are then calculated by 
solving the equations (39) and (40) simultaneously. 

In order to estimate the cluster physical parameters at r^tjo we 
use the definition of gas concentration parameter C500 to estimate 
'■500 > 

^•500 



C500 - 



(43) 



Wefixcsoo to the value given in Amaud et al. (2010) (C500 = 1.156). 
With knowledge of f-500 and all the four normalisation coefficients 
(Pei, ^ei, kjiTa, nd) we can calculate MT(r5oo), Mg(r5oo), /g(r5oo) and 
kETsoo. 



7 SIMULATED AMI SA DATA 

In generating simulated SZ skies and observing them with a 
model AMI SA, we have used the methods outUned in Hobson & 
Maisinger (2002) and Grainge et al. (2002). 

Generating a simulated cluster SZ signal using the isothermal 
j8-model requires the input parameters of z, T^, ne(0), r^, and j8; this 
set of parameters fully describes the Comptonisation y parameter. 
However, in order to verify the results of our analysis and to see if 
our methodology is capable of recovering the true values associated 
with the simulated cluster, it is instructive to estimate the cluster 
physical parameters using the three parameterisation discussed. We 
note that any parameterised model within the hierarchical structure 
formation of the universe for the ICM, including the isothermal 
jS-model, introduces constraints and biases in the inferred cluster 
parameters. Moreover, we now show that it is possible to get dif- 
ferent cluster physical parameters with the same set of input model 
parameters derived using the two different mass-temperature rela- 
tions described in Section 6.1.2. 

For example, if we consider parameterisation II we can use the 
HSE mass-temperature relation to calculate r2oo given in equation 
(29). Mg(r2oo) and MT(r2oo) are then calculated applying spherical 
geometry assumptions described in equations (24) and (25) respec- 
tively. We can also determine /g(r2oo)- 

However, if we consider parameterisation III we first calculate 
^t(''20o) using the virial mass-temperature relation given in equa- 
tion (33). r2oo and Mg(r2oo) are then estimated assuming spherical 
geometry for the cluster, (equations 25 and 24). A numerical exam- 
ple that leads to different results in cluster parameters is given in 
Table 5. To address this issue, we studied how /g(r2oo) varies as a 
function of and p while «e(0), Tg and z are fixed for both the HSE 
and virial M-T relations. Clearly, to obtain consistent results from 
both parameterisations one should select the values of rc and ji for 
which the corresponding gas mass fraction ratio is one. 
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Table 5. An example of input cluster parameters for the isothermal ,8-model 
that lead to inconsistent results using different parameterisations assuming 



h = 0.7. 


Input 


Assumed 


Derived 


Parameterisation 


Parameterisation 


parameter 


value 


parameter 


i,n 


III 


Xc 





''200 


1.56 Mpc 


1.71 Mpc 


yc 











rc 


200 kpc 




5.83X lO'^Mo 


7.67 X lO'^^Mo 


P 


0.95 










5keV 




5.91 X 10" Mg 


6.28 X 10" Mo 


ne(0) 


lO'^m-^ 








Z 


0.3 


fgirioo) 


0.102 


0.082 



Table 6. Cluster parameters for the isothermal yS-model assuming h = 0.7. 



Input 


Assumed 


Derived 


Parameterisation 


Parameter 


value 


Parameter 


I,II,III 







''200 


1.56 Mpc 


yc 





A^t('"200) 


5.83x101'' Mo 


rc 


155 kpc 


Mg('-200) 


6.36x10" Mo 


P 


0.79 


/g(''200) 


0.109 




5keV 


''500 


0.98 Mpc 


ne(0) 


lO'^m-s 


Mt(?-5oo) 


3.64 X 10''* Mo 


Z 


0.3 


Mg(r5oo) 


4.13 X 10"Mo 






/g(''500) 


0.11 



It should be noted that the same study may be carried out by 
investigating the variation of the ratio of gas mass fractions with 
either ne(0) or Tg while keeping rc and jS constant. However we find 
that the ratio is not sensitive to variation of these two parameters. 

Given the above, we decided to generate a simulated cluster 
with input parameters given in Table 6, which leads to consistent 
physical parameters for the cluster, at both r2oo and t^qq in both 
parameterisations. We assume that our cluster target is at declina- 
tion 6 = +40° observed for hour angles between - 4 and + 4 with 
2-s sampling for four days and 8 hours per day. We calculate the 
Comptonisation y parameter on a grid of 512 x 512 pixels with 
pixel size of 30". A realisation of the primary CMB is calculated 
using a power spectrum of primary anisotropics which was gen- 
erated for / < 8000 using CAMB (Lewis, Challinor & Lasenby 
2000), with a ACDM cosmology: Qm = 0.3, Ha = 0.7, erg = 0.8, 
Hd = VOkms^'Mpc"', Wq = -1 and Wa = 0. The CMB realisation 
is then co-added to the cluster in brightness temperature. It should 
be noted that in our simulation we did not include extragalactic ra- 
dio sources, or diffuse foreground emission from the galaxy as we 
have already addressed the efl^ects of the former in Feroz et al. 2009 
and the foreground galactic emission is unlikely to be a major con- 
taminant since our interferometric observations resolve out such 
large scale emission. The map is scaled by the primary beam ap- 
propriate to the measured value in the frequency channel and trans- 
formed into the Fourier plane. The resulting distribution is sam- 
pled at the required visibility points and thermal noise of 0.54 Jy 
per channel per baseline in one second which is appropriate to the 
measured sensitivity of the SA is added. Fig. 1 shows a map of the 
SZ temperature decrement of the first simulated cluster generated 
using the isothermal >8-model. 

To generate the second simulated cluster we use the GNFW 
pressure profile to calculate the Comptonisation y parameter. The 
input parameters for this model, (Pet. 'p), were selected to represent 
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Figure 1. Map of the simulated AMI SZ temperature decrement generated 
with the isothermal-jS model and parameters given in Table 6. Contours at 
... - 3, -2, 2, 3... times noise (cr = 79 (tJy beam"'), negative contours are 
dashed. The coordinates are J2000.0. 



Table 7. Cluster parameters for the entropy-GNFW model assuming h = 
0.7. 



Input 


Assimied 


Derived 


"Entropy"-GNFW 


parameter 


value 


parameter 


model 


Xc 





'"200 


1.56 Mpc 


yc 





*^t(''20o) 


5.83 X 10''' Mo 


'■p 


0.85 Mpc 


Mg(r2oo) 


6.36 X lO'^Mo 




37647.51 keVm-3 


/g(''200) 


0.109 


Z 


0.3 




0.98 Mpc 


rc 


155 kpc 




3.64 X lO'^^Mo 


a 


0.55 


Mg(r5oo) 


4.08 X 10" Mo 


Tim 


5keV 


./a(''5()()) 


0.11 








6.95 kcV 



a cluster with the same physical parameters at rioa and the same 
noise level as the first cluster. Although the parameters rc, a, z, and 
Tc{r2W)) do not contribute to the calculation of Comptonisation y 
parameter directly, their values were used to derive the parameters 
describing the GNFW pressure profile by following the steps de- 
scribed in section 6.2 to ensure that they represent the cluster with 
required physical parameters at r2oo. The cluster physical parame- 
ters at different radii will be difl'erent from the first cluster due to 
the different models describing the ICM and relaxing the assump- 
tion of isothermality in the second cluster. A summary of the cluster 
parameters is presented in Table 7. Fig. 2 shows a map of the SZ 
temperature decrement of the second simulated cluster, generated 
using the GNFW model. 



8 ANALYSIS AND RESULTS 

In this section we present the results of our analysis for all three 
parameterisations within the context of the isothermal yS-model and 
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Figure 2. Map of the simulated AMI SZ temperature decrement gener- 
ated with GNFW model and parameters given in Table 7. Contours at 
... - 3,-2,2,3... times noise (cr = 81.5 |J.Jy beam"'), negative contours 
are dashed. The coordinates are J2000.0. 



the entropy-GNFW pressure model. In each case we first study our 
methodology in the absence of data. This can be carried out by 
setting the likelihood to a constant value and hence the algorithm 
explores the prior space. This analysis is crucial for understanding 
the underlying biases and constraints imposed by the priors and the 
model assumptions. Along with the analysis done using the simu- 
lated AMI data, this approach reveals the constraints that measure- 
ments of the SZ signal place on the cluster physical parameters and 
the robustness of the assumptions made. It should be noted that in 
all the plots of probability distributions, we explicitly include the 
dimensionless Hubble parameter h = /fo/(100kms"'Mpc"') with 
h set to 1.0. 



8.1 Analysis using isothermal /J-model-Parameterisation I 

Figs. 3 and 4 represent the results of a priors-only analysis 
showing the sampling and derived parameters respectively. 1-D 
marginalised posterior distributions of sampling parameters in Fig. 
3 show that we were able to recover the assumed prior probability 
distributions for cluster position and the gas mass. However, this 
parameterisation clearly prefers higher temperature and yS and the 
probability distribution for falls as we go towards higher r^. This 
feature in particular creates a void region in the 2-D marginalised 
probability distributions of yS - and Tg - at higher which im- 
plies that low mass clusters are unlikely to have high and low j}. 
This effect is a direct result of imposing the constraint that rsm > 0. 
Moreover, as may be seen from Fig. 4, this choice of priors drives 
the posterior probability distributions of both the gas mass and the 
gas mass fraction towards low values. 

Figs. 5 and 6 show the results of the analysis of the simulated 
SZ cluster data. The vertical lines show the true values of the pa- 
rameters. Table 8 also summarises the mean, the dispersion and the 
maximum likelihood of each parameter. 

In Fig. 5, we notice the strong degeneracy between and /? 




Figure 3. Marginal distiibutions for the sampling parameters with no data 
for isothermal j8-model -parameterisation I. 
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Figure 4. Marginal distributions for the derived cluster physical parameters 
with no data for isothermal /J-model -parameteiisation I. 



(Grego et al. 2001). However, it is apparent that neither p nor T„ 
is well-constrained using this parameterisation. Also, higher values 
than the true input parameters are preferred for both parameters. 
This effect leads to two results: firstly it yields a higher estimate for 
r2oo and so equation (28) overestimates the total mass; secondly, 
since for this parameterisation there is a negative degeneracy be- 
tween gas mass and temperature, the high temperature therefore 
leads the marginalised posterior distribution for gas mass peaking 
towards the lower end of the distribution although the recovered 
mean value of Mg(r2oo) is within Icr from its corresponding input 
value for the simulated cluster. As a result of these two effects, the 
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Figure 5. Marginal distributions of tlie sampling parameters for simulated 
data for isothermal /3-model -parameterisation I. 
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Figure 6. Marginal distributions of the derived cluster physical parameters 
with simulated data for isothermal /3-model - parameterisation I. 

gas mass fraction is driven even further to the lower end of the 
allowed range. There is also a degeneracy between the two free 
parameters of jS and Tg; this degeneracy again originates from de- 
pendency of r2oo on both parameters as given in equation (29). 

8.2 Analysis using isothermal yS-model-Parameterisation II 

Figs. 7 and 8 show the results from prior-only analysis for parame- 
terisation II. We recover the assumed prior probability distributions 
for cluster position, f3, total mass and gas mass fraction. There is 
a similar trend in the 1-D posterior probability distribution of to 



Table 8. Simulated cluster parameters (mean, standard deviation and Max- 
imum likelihood) estimated using isothermal yS-model-parameterisation I 



assuming h 


= 0.7. 




1 arameter 




M 




-2.1 ± 15.3" 


-0.6" 


yo 


6.3 ± 13.9" 


5.5" 




391.26 ± 214.98 kpc 


138.77 kpc 


p 


1.7 ±0.5 


0.98 


Mg(r2oo) 


(3.41 ± 3.16) X lO'^Mo 


2.86 X lO''' Mo 




10.61 ±5.28 keV 


10.15 keV 


AfT(r5oo) 


(3.91 ± 3.41) X 10'^ Mo 


1.49 X lO'^Mo 


^500 


1.96 ± 0.69 Mpc 


1.57 Mpc 


Mg(r5oo) 


(2.71 ± 1.42) X 10'' Mo 


2.31 X lO'^Mo 


/g('-500) 


0.15 ±2.7 


0.014 


Mi(r2m) 


(6.43 ± 5.43) X 10'5 Mo 


2.37 X 10'^ Mo 




3.14 ± 1.0 Mpc 


2.49 Mpc 


fginoo) 


(0.14 ±3.4) 


0.012 



Table 9. Simulated cluster parameters (mean, standard deviation and Max- 
imum likelihood) estimated using isothermal /?-model-parameteiisation II 



assuming h 


= 0.7. 




Parameter 


p ± cr 






-2.6 ± 15.7" 


5.5" 


yc 


6.4 ± 14.5" 


5.6" 




410.37 ±237.24 kpc 


135.43 kpc 


P 


1.7 ±0.5 


0.8 


Mi(r2m) 


(6.8 ± 2.1) X 10i'*Mo 


5.0 X 10''' Mo 


/g(''200) 


0.12 ±0.03 


0.11 


MT(r5oo) 


(3.5 ± 8.81) X 10" Mo 


3.13 X W*Mq 


'■500 


0.96 ±0.08 Mpc 


0.93 Mpc 


Mg(r5oo) 


(6.2 ± 1.6) X lO'^Mo 


3.9 X lO'-'' Mo 


/g(''500) 


0.18 ±0.05 


0.12 


''200 


1.59 ± 1.57 Mpc 


1.47 Mpc 


A?e(''200) 


(7.76 ± 2.08) X 10'')Mo 


5.35 X lO'^Mo 




3.0 ± 1.2keV 


4.3 keV 



that mentioned in the parameterisation I, which leads to a void re- 
gion in the 2-D marginalised posterior distribution of M-^irooQ) — 
for the same reason as discussed for the parameterisation I. How- 
ever, parameterisation II prefers a lower temperature which arises 
from the fact that HSE mass-temperature relation used in this pa- 
rameterisation (equation 32) is inversely proportional to /?. 

Figs. 9 and 10 show the results of the analysis using simu- 
lated SZ cluster data, with vertical lines representing the true pa- 
rameter values. Table 9 also summarises the mean, the dispersion 
and the maximum likelihood values of each cluster parameter es- 
timated using parameterisation II. A tight degeneracy between /"c 
and p is noticeable in the corresponding 2-D marginalised proba- 
bility distribution. /3 on the other hand is not well constrained and 
moves towards higher values which results in the probability distri- 
bution of temperature being driven to lower values again because of 
the 1 IP relationship in equation (32). However, this parameterisa- 
tion along with the simulated SZ data reliably constrains Mt(?'2oo), 
Wa(''20o) and /g(r2oo). Comparing the 1-D marginalised posterior 
distributions of gas mass fractions at two overdensity radii rsoo and 
rioo also reveals that we cannot constrain the radial behaviour of 
the gas mass fraction using this parameterisation, as /g(r5oo) ex- 
hibits too wide a probability distribution. For /g(r2oo), we seem to 
have recovered the input prior distribution. 
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Figure 7. Marginal distributions of tlie sampling parameters with no data 
for isothermal j8-model- parameterisation II. 
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Figure 8. Marginal distributions of the derived cluster physical parameters 
with no data for isothermal /3-model-parameterisation II. 
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Figure 9. Marginal distributions of the sampling parameters with simulated 
data for isothermal /3-model- parameterisation II. 
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Figure 10. Marginal distributions of the derived cluster physical parameters 
with simulated data for isothermal /3-model-parameterisation II. 



8.3 Analysis using isothermal yS-model-Parameterisation III 

The results of the analysis with no data are plotted in Figs. 1 1 and 
12. It is evident that, while the assumed prior probability distribu- 
tions for the cluster position, total mass and gas mass fraction are 
recovered, the two sampling parameters and /? show the same be- 
haviours as discussed for the other two parameterisations. We also 
see a trend towards lower values in the 1-D posterior probability 
distribution of temperature. However this behaviour is due to the 
direct relationship between the total mass and the temperature in 
this parameterisation and the specific prior distribution we have as- 



sumed for the total mass which clearly has a higher probability at 
the lower masses. 

Figs. 13 and 14 represent the marginalised posterior distribu- 
tions from the analysis of simulated SZ cluster data for sampling 
and derived parameters respectively while in Table 10 we present 
the mean, the dispersion and the maximum likelihood of each clus- 
ter parameter estimated using parameterisation III. The strong de- 
generacy between and yS is quite apparent in this parameterisa- 
tion, while p is poorly constrained and biased towards higher val- 
ues. We note that, since the SZ analysis constrains cluster total mass 
internal to the radius ^200 and we use the virial M-T relation (equa- 
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Figure 11. Marginal distributions of the sampling parameters with no data 
for isothermal /?-model- parameterisation III. 



Figure 13. Marginal distributions of the sampling parameters with simu- 
lated data for isothermal yS-model- parameterisation III. 
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Figure 12. Marginal distributions of the derived cluster physical parameters 
with no data for isothermal /3-model- parameterisation III. 



Figure 14. Marginal distributions of the derived cluster physical parameters 
with simulated data for isothermal /3-model- parameterisation III. 



tion (33)) to derive cluster average temperature within this radius, 
the result of temperature estimation is less biased and more reli- 
able than the parameterisations I and II in recovering the tempera- 
ture true value. We have used this parameterisation in our follow-up 
analysis of the real data where we studied a joint weak gravitational 
lensing and SZ analysis of six clusters (AMI Consortium: Hurley- 
Walker et al. 2011) and high and moderate X-ray luminosity sam- 
ple of LoCuSS clusters (AMI Consortium: Rodriguez-Gonzalvez 
et al. 201 1; AMI Consortium: Shimwell et al. 201 1). 



8.4 Analysis using entropy-GNFW pressure model 

Similar to the isothermal y8-model we first studied our methodology 
for the "entropy"-GNFW pressure model with no data. The results 
are represented in Figs. 15 and 16. This analysis again helps us 
understand which parameters are constrained by SZ measurement 
as well as to check the algorithm in retrieving the prior probabil- 
ity distributions. From both 1-D and 2-D marginalised probability 
distributions it is clear that we are able to recover the input pri- 
ors probability distributions and the probability distributions of the 
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Table 10. Simulated cluster parameters estimated (mean, standard deviation 
and Maximum likelihood) using isothermal /3-model-parameterisation III 



assuming h 


= 0.7. 




parameter 


H±cr 


/J 




-3.0 ± 15.6" 


-5.0" 




6.4 ± 14.4" 


5.7" 


rc 


395.11 ± 226.21 kpc 


142. 14 kpc 


P 


1.7 ±0.5 


1.1 


MT(r2oo) 


(4.68 ± 1.56) X 10'* Mq 


4.46 X 10"* Mo 




0.11 ±0.03 


0.1 




(7.35 ± 3.0) X 10"* Mq 


4.58X lO'^^Mo 




1.21 ±0.17 Mpc 


1.06 Mpc 


Mg(r5oo) 


(4.50 ± 1.04) X IO'^^'Mo 


4.00 X lO'"* Mo 


/g(''500) 


0.07 ± 0.03 


0.09 


'"200 


1.43 ± 1.50 Mpc 


1.42 Mpc 


Mg(r2oo) 


(5.14 ± 1.6) X 10" Mo 


4.49 X 10'^^ Mo 




4.3±0.9keV 


4.2 keV 



derived parameters are according to their corresponding functional 
dependencies on the sampling parameters. 

Figs. 17, 18 and Table 11 show the results of our analysis for 
"entropy" -GNFW pressure profile using parameterisation II while 
Figs. 19, 20 and Table 12 show the results of the same analysis 
using parameterisation III. We note that in both analyses and a 
the parameters that define the shape of the entropy profile are not 
constrained while the scaling radius, r^, which defines the GNFW 
pressure profile is completely constrained. As a result we notice 
similar constraints in the estimation of rjoo in both parameterisa- 
tions since we assume a fixed cjoo- We also note the degeneracies 
between Mx(r2oo)-''c and MT(r2oo)-Q' which are because of the de- 
pendency of Pei on these two free parameters. On the other hand the 
^t(''20o)- '■p degeneracy seen in Figs. 17 and 19 is due to the intrin- 
sic degeneracy that exists between the cluster size and the volume 
integrated Comptonisation parameter (Y^z-f^ degeneracy) in the SZ 
measurements (Planck Collaboration 201 Id). Moreover, compar- 
ing rg(r5oo) and rg(r2oo) (Table 1 1 and 12) confirms a radial decline 
in the ICM temperature distribution as expected. 

Overall, both parameterisations could constrain the cluster 
physical parameters, however, analysis using parameterisation III 
leads to a tighter constrain on both r„(rjoo) and rg(r2oo)- The results 
of Parameterisation III once more show that this parameterisation 
can reliably be used in the analysis of clusters of galaxies as it is 
less model dependent and produces unbiased results in particular 
when the assumption of hydrostatic equilibrium breaks, in young 
or disturbed clusters (parameterisation II). 




,5 1.5 2.5 2 4 0.050.10.15 1 2 3 



Figure 15. Marginal distributions of the sampling parameters with no data 
for "entropy"- GNFW pressure model. 






9 DISCUSSION AND CONCLUSIONS 

We have studied two parameterised models, the traditional isother- 
mal /3-model and the "entropy"-GNFW pressure model, to analyse 
the SZ effect from galaxy clusters and extract their physical pa- 
rameters using AMI SA simulated data. In our analysis we have 
described the current assumptions made on the dynamical state of 
the ICM including spherical geometry, hydrostatic equilibrium and 
the virial mass-temperature relation. In particular we have shown 
how different parameterisations which relate the thermodynamical 
quantities describing the ICM to the cluster global properties via 
these assumptions lead to biases on the cluster physical parameters 
within a particular model. 

In this context, we first generated a simulated cluster using 



Figure 16. Marginal distributions of the derived cluster physical parameters 
with no data for "entropy"-GNFW pressure model. 



the isothermal y6-model observed with the AMI SA and used these 
simulated data to study three different parameterisations in deriv- 
ing the cluster physical parameters. We showed that in generating 
AMI simulated data, it is extremely important to select the model 
parameters describing the SZ signal in a way that leads to the con- 
sistent cluster parameter inferences upon using the three different 
parameterisation methods. 

We found that each parameterisation introduces different con- 
straints and biases in the posterior probability distribution of the 
inferred cluster parameters which arise from the way we imple- 
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Table 11. Simulated cluster parameters estimated (mean, standai'd deviation 
and Maximum likelihood) using "entropy"-GNFW pressure model- param- 



eterisation II assuming h = 0.7. 


Paramctc r 


/j_o- 






-8.9 ± 10.04" 


-9.3" 




6.3 ± 9.5" 


7.8" 




770.43 ± 3o2.29kpc 


1355.86 kpc 


a 


1.2 ±0.67 


1.2 


MT(r2oo) 


(5.86 X 10''*±3.43x 10")Mo 


3.3 X 10" Mo 


/g('-200) 


0.11 ±0.02 


0.13 


''p 


1.03 ± 0.33 Mpc 


0.87 Mpc 




(8.86 ± 11.57)x 10'''Mo 


4.0 X 10" Mo 


'■500 


1.18 ±0.39 Mpc 


1.01 Mpc 




(5.71 ±4.08) X 10'' Mo 


3.06 X lO'^Mo 


fgirsoo) 


0.1 ±0.09 


0.07 


Tg(rsoo) 


7.5 ± 1.8 keV 


6.5 keV 


''200 


1.57 ±0.28 Mpc 


1.3 Mpc 




(6.53 ± 3.67) X IO'^Mg 


4.14 X 10" Mo 


Tg{r2oo) 


7.4±2.6keV 


5.02 keV 



Figure 17. Marginal distributions of the sampling parameters with simu- 
lated data for "entropy"-GNFW pressure model using parameterisation 11. 
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Figure 18. Marginal distributions of the derived cluster physical parameters 
with simulated data for "entropy"-GNFW pressure model using parameter- 
isation II. 



ment assumptions about the cluster structure and its composition. 
The biases in the posterior probability distributions of the cluster 
parameters are more pronounced in parameterisations I and II, as 
the results depend strongly on the relatively unconstrained clus- 
ter model shape parameters: and yS. However, the biases intro- 
duced by the choice of priors are even worse in parameterisation 
I, in which the gas temperature is assumed to be an independent 
free parameter. This, along with the assumption of isothermality, 
causes the priors to dominate in extracting the cluster physical pa- 
rameters regardless the type of prior chosen for the gas tempera- 



Table 12. Simulated cluster parameters estimated (mean, standard devi- 
ation and Maximum likelihood) using entropy-GNFW pressure model- 
parameterisation 111 assuming h = 0.7. 

Parameter ii±(t u 



Xq 


-8.8 ±9.8" 


-10.4" 


yc 


6.6 ± 9.5" 


7.2" 


ro 


798.7 ±376.8 kpc 


1280.8 kpc 


a 


1.3 ±0.69 


1.3 




(8.0 X 10"'±5.6x 10''')Mo 


4.5 X IO'-'Mo 


fgirioo) 


0.11 ±0.03 


0.13 


rp 


1.0 ±0.3 Mpc 


0.86 Mpc 




(6.6 ± 1.14) X 10"* Mo 


3.4 X 10'" Mo 


'■500 


1.14 ±0.43 Mpc 


1.0 Mpc 


A'g('-500) 


(7.35 ± 5.92) X lO'^Mo 


3.67 X lO'^Mo 


/gCf-soo) 


0.14 ±0.14 


0.11 




6.0 ± 1.4keV 


5.5 keV 




1.7 ±0.28 Mpc 


1.4 Mpc 




(9.18 ± 5.71) X 10" Mo 


5.91 X lO'^Mo 


Tg{r2Qo) 


5.98 ± 2.43 keV 


4.2 keV 



ture (AMI Consortium: Rodrfguez-Gonzalvez et al. 2011 and AMI 
Consortium: Zwart et al. 2010). The cluster physical parameters 
estimated using parameterisation I depend strongly on the model 
parameters. Although it can constrain the cluster position and its 
Mj(r2oo), it fails to recover the true input values of most of the sim- 
ulated cluster properties. For example the inferred values for mass 
and temperature at r2oo are Mj{r2oo) = (6.43 ± 5.43) x lO'"" Mo 
and Tg(r2Qo) = (10.61 ± 5.28) keV whereas the corresponding in- 
put values of simulated cluster are: MT(r2oo) = 5.83 x lO''* Mo 
and T^(r2oo) = 5keV. In terms of the application to the real data, 
we have noticed similar biases in the results of our analysis of 7 
clusters using this parameterisation (AMI Consortium: Zwart et al. 
2010). In order to improve our analysis methodology in parame- 
terisations II and III, the correlation between the cluster total mass 
and its gas temperature is taken into account. In parameterisation 
II we relate MT(r2oo) and rg(r2oo) using the hydrostatic equilibrium 
whereas in parameterisation III we use virial mass-temperature re- 
lationship. It should be noted that the derived rg(r2oo) in parameter- 
isation II is the gas temperature at the overdensity radius r2oo which 
is then assumed to be constant throughout the cluster. In parameter- 



© 201 1 RAS, MNRAS 000, 1-19 



16 



AMI Consortium: Olamaie et al. 




Table 13. The results of 100 CMB realisations for the thi'ee parameterisa- 
tions assuming h=0.7. 
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Figure 19. Marginal distributions of the sampling parameters with simu- 
lated data for "entropy"-GNFW pressure model using parameterisation III. 
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Figure 20. Marginal distributions of the derived cluster physical parameters 
with simulated data for "entropy"-GNFW pressure model using parameter- 
isation III. 



isation III, however, rg(r2oo) is the mean gas temperature internal to 
radius r2oo and is assumed to be constant. We notice that analysing 
the same simulated data set using parameterisation II can constrain 
the 1-D posterior distribution of the cluster physical parameters bet- 
ter than parameterisation I such that MT(r2oo) = (6.8±2.1)xl0''' Mq 
and T„{r2m) = (3.0 ± 1.2) keV. Since parameterisation II uses the 
full parametric hydrostatic equilibrium, the temperature estimate 
depends on and f3 and is therefore biased low. These results were 
also confirmed in our analysis of the bullet like cluster A2146 (AMI 
Consortium: Rodriguez-Gonz alvezet al. 201 1). Relating the cluster 



parameterisation 


MT(r2oo) Mo 


7'g(r2oo)keV 


I 


(6.18 ± 5.23) X 10'5 


11.18±5.16 


n 


(8.067 ± 2.61) X lO'"* 


3.94 ± 1.67 


III 


(5.94 ± 2.26) X 10"* 


4.97 ± 1.21 



total mass and its temperature via virial theorem in parameterisa- 
tion III leads to less bias in cluster physical parameters compared 
to the other two parameterisations as it is less model dependent: 
MT(r2oo) = (4.68 ± 1.56) x 10" Mq and rg(f-2oo) = (4.3 ± 0.9) keV. 

A detailed comparison between our different parameterisa- 
tions both using simulated data and on the bullet like cluster A2146 
(AMI Consortium: Rodriguez-Gonzalvez et al. 2011) found that 
parameterisation III can give more reliable results for cluster phys- 
ical properties as it is less dependent on model parameters. Param- 
eterisation II also gives convincing estimates for the cluster total 
mass and its gas content although its temperature estimate is poorly 
justified, as it depends strongly on the model parameters. Moreover, 
young or disturbed clusters are unlikely to be well-described by 
hydrostatic equilibrium. We therefore used parameterisation III as 
our adopted analysis methodology in our follow-up studies of the 
real clusters including the joint SZ and weak lensing analysis of 
six clusters (AMI Consortium: Hurley- Walker et al. 2011) and the 
analysis of LoCuss cluster sample (AMI Consortium: Rodriguez- 
Gonz alvez et al. 2011; AMI Consortium: Shimwell et al. 201 1). 

In order to make sure that our results are not biased by one 
realisation of primordial CMB, we have studied 100 CMB realisa- 
tions for the three parameterisations. The 1-D marginalised poste- 
rior probability distributions of MT(r2oo) and Tg(r20Q) are shown in 
Figs. 21, 22 and 23 for each parameterisation. The solid blue line 
represents the true value corresponding to the simulated cluster and 
the dashed red line shows the mean value of the distributions. Table 
13 also presents the numerical results of this analysis. Comparing 
the ID posterior distributions along with the mean values of the 
MT(r2oo) estimates in the three parameterisations for these 100 real- 
isations show that parameterisation I can hardly constrain the sim- 
ulated cluster properties and recover the input true values. Param- 
eterisation II can constrain the cluster total mass, however, the gas 
temperature estimate is biased low as it depends on unconstrained 
model shape parameters. On the other hand, parameterisation III 
can indeed constrain both cluster mass and its gas temperature and 
the results are unbiased. 

In order to remove the assumption of isothermality which is of 
course a poor assumption both within the cluster inner region and 
at the large radii and to improve our analysis model for the clus- 
ter ICM which can be fitted accurately throughout the cluster, we 
also studied the SZ effect using "entropy"-GNFW pressure model. 
This model assumes a 3-D /J-model like radial profile describing 
the entropy in the ICM as well as the GNFW profile for the plasma 
pressure. This choice is reasonable as the entropy is a conserved 
quantity and describes the structure of the ICM while the pressure 
is related to the dark matter component of the cluster. Moreover, 
among all the thermodynamical quantities describing the ICM, en- 
tropy and pressure show more self- similar distribution in the out- 
skirts of the cluster. The combination of these two profiles then 
allows us to relate the SZ observable properties to the cluster phys- 
ical parameters such as its total mass. This model also allows the 
electron pressure and its number density profiles to have different 
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Figure 21. 100 realisations of 1-D marginalised posterior probability 
distributions of -Mt(''20o) and 7'g(''20o) using isothermal /beta-model- 
parameterisation I. The solid blue line represents the true value correspond- 
ing to the simulated cluster and the dashed red line shows the mean value 
of the distributions. 



Figure 22. 100 realisations of 1-D marginalised posterior probability dis- 
tributions of MT(r2oo) and 7'g(''20o) using isothermal /?-model- parameter- 
isation II. The solid blue line represents the true value corresponding to 
the simulated cluster and the dashed red line shows the mean value of the 
distributions. 



distributions leading to a 3-D radial temperature profile. In this con- 
text we simulated a second cluster using an entropy-GNFW pres- 
sure profile with the same physical parameters and thermal noise as 
the first cluster at ^200- 

We then analysed the second simulated cluster using 
"entropy"-GNFW pressure model with different parameterisations. 
In this model temperature is no longer isothermal so that we can 
not use parameterisation I where a single temperature is assumed 
as an independent input parameter. The results of our analysis using 
parameterisation II and III show that while the characteristic scal- 
ing radius describing the GNFW pressure profile is constrained, 
the shape parameters defining the entropy profile remain uncon- 
strained. Moreover, all the cluster physical parameters lie within 
Icr errorbars from the corresponding true values of the simulated 
cluster in the two parameterisations. However, parameterisation III 
provides tighter constrains in 1-D marginalised posterior distribu- 
tion of the temperature and the overall results are less model de- 
pendent so that it can be reliably used in the analysis of galaxy 
clusters in particular when the assumption of hydrostatic equilib- 
rium breaks (e.g. in disturbed clusters and clusters that are going 
through merging). 

We conclude that using the "entropy"-GNFW pressure model 
overcomes the limitations of the isothermal yS-model in fitting clus- 
ter parameters over a broad radial extent. However, AMI simulated 
data do not strongly prefer one model over the other. We investi- 
gated this conclusion further by fitting both GNFW pressure profile 
and isothermal /3-model to a simulated cluster with 6500 = 2.5' and 
Kjoo = 2.5 X 10"'(arcmin)-. The result is shown in Fig. 24 with blue 
dashed line representing the fit using the isothermal j8-model and 
the red representing the fit using GNFW pressure profile. However, 
we aim to compare these two models in our future studies using the 
real data. 
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Figure 23. 100 realisations of 1-D marginalised posterior probability dis- 
tributions of MT(r2oo) and 7'5(r2oo) using isothermal /J-model- parameter- 
isation III. The solid blue line represents the true value corresponding to 
the simulated cluster and the dashed red line shows the mean value of the 
distributions. 
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Figure 24. The SZ flux amplitude versus AMI SA observing baseline for a 
cluster with 6500 = 2.5' and {500 = 2.5 X 10"'(arcmin)~. Blue dashed line 
represents the fit using the isothermal jS-model and the red represents the fit 
using GNFW pressure profile. 
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